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Abstract 

In an attempt to uncover any underlying structure in the neutrino mass matrix, 
we discuss the possibility that the ratios of elements of its Major ana mass matrix 
are equal. We call this "strong scaling Ansatz" for neutrino masses and study its 
phenomenological implications. Of three possible independent scale invariant possi- 
bilities, only one is allowed by current data, predicting in a novel way the vanishing 
of t/e3 and an inverted hierarchy with the lightest neutrino having zero mass. The 
Ansatz has the additional virtue that it is not affected by renormalization running. 
We also discuss explicit models in which the scaling Ansatz is realized. 



*email: rmoliapat@pliysics.umd.edu 
t email: werner_rode j ohannOph . turn . de 



1 Introduction 



Understanding the new physics behind the neutrino observations is one of the main chal- 
lenges for theory [1] today. While it is widely expected that the extreme smallness of the 
neutrino masses is most likely related to the seesaw mechanism [2], i.e, the existence of 
a B — L symmetry broken at very high scale, there is no common consensus on how to 
understand the detailed mixing pattern. Observation implies a particular structure in the 
neutrino mass matrix 
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M u = UMf ag U T , (1) 



whose diagonalization in the charged lepton mass basis gives the PMNS matrix U. There 
are several distinct approaches: 

(i) symmetry approach, where one postulates family symmetries to constrain the mass 
matrix elements and studies their consequence; 

(ii) anarchy approach, where all elements of the neutrino mass matrix M. v are allowed to 
vary at random and the most probable values are then confronted with experiments 
and 

(iii) top-down GUT approach, where one considers a quark-lepton unified grand unified 
framework such as 5*0(10) and obtains predictions for neutrino parameters. 

In this paper, we consider a different approach which is located somewhere between the 
anarchy and the symmetry approaches. We consider the Ansatz that the elements of the 
neutrino mass matrix, in the basis in which the charged lepton mass matrix is diagonal, 
obey scaling laws. According to this law, ratios of certain elements of M.ij are same. We 
will call this the "strong scaling Ansatz" (SSA) on neutrino masses. An important property 
of the SSA is that it is not affected by renormalization group extrapolation as we run it 
down from the seesaw scale to the weak scale. This property is not shared by many other 
Ansatze used in literature such as fi-r or other symmetries, texture zeros, etc. 
The first question to ask now is which matrix elements should one impose the scaling 
Ansatz to? If we consider scaling for a given matrix element as m a p ~ m r Xa+x i 3 , where 
mo is flavor independent, then this implies three kinds of strong scaling Ansatze of the 
form that — — is independent of the flavor oti. In terms of matrix elements, this means 

(A) : ^ = ^tH = VhE = c . (2) 

Tfl eT 1TI[_it Tfl TT 

and cyclic variations of these relations, i.e., 

(B) ' trifle ITT're 

TTlvT tTtlin- TYItt 
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and 

(UJ : = = = C . (4J 

TTiefi TTlfin i^Wfj, 

These are the only three independent strong scaling Ansatze for the 3x3 case. In the 
next Section [2j we will consider the phenomenological implications of these strong scaling 
constraints, showing in particular that only case (A) is allowed by current data. 
Clearly, we do not know the fundamental origin of the SSA, but if one looks at the hierarchy 
of fermion masses, such flavor dependent power laws for masses could easily exist and we 
only need to make some assumptions about the coefficients. Alternatively, SSA could arise 
from seesaw physics, i.e., from the structure of the Dirac and heavy Majorana neutrino 
mass matrices. In any case, we will discuss theoretical scenarios where SSA can emerge in 
Section [3j 



2 Phenomenology of the Strong Scaling Ansatz 

In order to explore the implications of SSA, we start by writing down the usual form of 
the PMNS matrix: 
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where P = diag(l, e ia , e l ^ + ^) contains the Majorana phases. Let us discuss the phe- 
nomenology of SSA, starting with case (A) defined in Eq. (j2J): 
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(6) 



To discuss the phenomenology of the model, we first note that the mass matrix in Eq 
has rank 2 and therefore predicts one vanishing neutrino mass. The eigenvector of Eq 
corresponding to the zero eigenvalue is given by (0, —1 



i + \c\ 2 ,c/^TTW) T - 



6) 

There- 



fore, the strong scaling condition is only compatible with an inverted mass hierarchy, 
U e 3 = and hence no CP violation in oscillation experiments. Thus, SSA provides an- 
other way to understand the vanishing of U e 3 and one can use small deviations from strong 
scaling as a way to understand any possible non-zero U e 3 value, see Section 12.31 for more. 
Atmospheric neutrino mixing is described by tan 2 #23 = | l/c| 2 , and data requires that |c| is 
close, but not necessarily equal, to 1. Let us assume for the moment that the parameters 
are all real. The non-zero masses read 
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The solar neutrino mixing angle is given by 
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and the effective mass governing neutrinoless double beta decay is m \ A\. To have the 
ratio AttIq/ Am 2 A = (m 2 — m^/mj ~ 1/25, typically D + c 2 (A + D) < w, or A ~ —2D, 
must hold. This is a tuning of parameters shared by many other Ansatze for neutrino mass 
models. In this case (with c = +1 one needs negative A) sin 2 6* 12 ~ | (1 + A/^A 2 + 2 £> 2 ) 
and therefore i? 2 > A 2 in order to have sin 2 #12 < |. This pattern of .M^ corresponds 
to the approximate conservation of the flavor symmetry L e — — L T [3], and one has 
(m) ~ a/ Am} cos 26>i2 = \/AmJ (1 - 2 sin 2 #12). There is also a second, somewhat 
tuned possibility for the parameters A,B,D, namely A 2 ,/} 2 > B 2 , which corresponds 
to (m) ~ a/ Am 2 A . 

One might compare SSA matrix to the case of exact /i-r symmetry [3], which also predicts 
Z/ e 3 = and in addition cos 2023 = 0. A mass matrix obeying /i-r symmetry reads 
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(9) 



Hence, /i-r symmetry with the additional constraint E = D reproduces SSA when in 
addition |c| = 1 holds. Breaking of /x-r symmetry will in general lead to non-zero U e ^ 
and cos 26*23, together with a correlation of these parameters which depends on the way 
the symmetry is broken [5]. SSA in turn has in general non-zero cos 26*23 and breaking of 
SSA (treated in Section 12.31) will generate non-zero ?7 e 3 not directly linked to the deviation 
from zero cos 26*23. Thus a key test that distinguishes this model from an approximate /i-r 
symmetric model is that here we can have departure from the maximal atmospheric mixing 
angle even though U e $ = 0, whereas in the case of generic breaking of /i-r symmetry, there 
is a strong correlation between 6*23 — f and U e3 . 

Another comparison should be done with the flavor symmetry L e — — L T , which is the 
usual symmetry to enforce an inverted hierarchy [3]. The mass matrix reads 
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(10) 



and predicts vanishing Am 2 & . In analogy to SSA, m 3 and U e ^ are zero and atmospheric 
neutrino mixing is non-maximal: tan 2 #23 = |A/5| 2 . However, solar neutrino mixing is 
maximal, and one requires large and tuned breaking of the symmetry in order to obtain a 
sufficiently non- maximal #12 [6]. The breaking parameters have to have at least 40% the 
size of the parameters allowed by the symmetry. In the case of SSA one has naturally large 
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but non- maximal solar neutrino mixing, as given in Eq. (jSj), and therefore does not face 
the problems of L e — — L T . 

It is interesting to elaborate on renormalization aspects of the strong scaling condition. 
Presumably, SSA will be a property introduced at the scale Mx ~ Mqut — 2 ■ 10 16 GeV 
of neutrino mass generation. Its predictions have to be compared to the measurements 
performed at low scale mz- It is well-known that the effect of running from high to low 
scale can be taken into account by multiplying the af3 entry of M. v with (1 + 5 a )(l + 6/3), 
where 5 a = C 1( ™i v i In ^j, with m a being the charged lepton mass [TJ. The parameter C 
is given by 3/2 in the SM and by -(1 + tan 2 /3) in the MSSM. Obviously, with 5 T > 6 £;fl , 
this leads in the strong scaling condition (T2J) only to a rescaling of c with (1 + 5 T ), i.e., 
the three ratios of the mass matrix elements still have the same value, which is however 
slightly changed. The predictions rn.3 = and U e3 = are therefore invariant under the 
renormalization running. The other predictions of SSA are affected: for instance, if SSA 
at high scale imprints the ratio c and therefore predicts tan 2 823 = 1/c 2 , then at low scale 
we have tan 2 823 = (1 — 2S T )/ c 2 . Note that \l-t symmetry and L e — L^ — L T are in general 
affected by renormalization running, hence one generically expects U e3 7^ for such models. 



2.1 CP violating case 

The fact that U e3 = means that there is no CP violation in oscillation experiments. This 
can also be seen by working with invariants: a useful measure of CP violation in neutrino 
oscillation is given by [8] 

J CP = I sin 28 12 sin 28 23 sin 28 13 cos 8 13 sin 5 = Im {^12 ^23^31} 

8 Am 21 Am 31 Am 32 (11) 

where h = M u Ml . 

From Eq. (J6j) it follows that ^13 = hyi/c* and /i23 = 0/133. Since h is hermitian, one has 

Im{h 12 h 23 h 31 } = Im {h 12 h 23 h* l3 } = Im {c* h 13 ch 33 h* 13 } Im { |c| 2 \h 13 \ 2 ti3 3 } = . (12) 

From hi3 = h\ 2 /c* and ^23 = ch 3 3 it also follows that ^12/^13 — ^32/^33 = 0. Using the 
parametrization of the PMNS matrix from Eq. (jSJ), this condition leads for an inverted 
mass ordering to [9] 

^ 1 |Am|m 2 sin26>i2 cot6> 23 

Am\ (ml + Am 2 A ) + Am| (Ami - m l sin2 ^12) 

Since we know that 7713 = 0, it automatically follows again that \U e3 \ = (and vice versa). 
Suppose now that all parameters in Eq. ([6]) are complex, i.e., A — > Ae 1 ^ 1 , B — > Be*, 
D — > D e 1 ^ 3 and c — > ce*^ 4 , where the redefined A, B, D, c are real. It is easy to show that 
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one can rephase the lepton fields in a way such that only the e/i and er entries are left 
with a phase = 2 — 0i/2 — 03/2: 
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(14) 
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There are thus five relevant parameters A, B, D,c,<p for the observables #12, 623, Am 
Am\ and (m). The Ansatz can therefore completely be reconstructed, in contrast to 
many other Ansatze. A 23 rotation with tan 823 = 1/c gives 
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(15) 



We note here that the SSA prediction of 7B3 = indicates that one of the Majorana phases, 
in this case (3, is absent. The only observable phase is therefore the Majorana phase a. 
With A, B, D typically of the same order of magnitude, it is rather cumbersome to rephase 
this simple looking matrix in order to identify the sole surviving low energy phase a and 
also the two non-vanishing masses. Again, invariants are helpful: the in general three low 
energy phases a, (3,5 correspond to three independent invariants, which can be chosen to 
be [TO] 

I af3 = hn.{m aa mi3pm* a pm^ a } , where (a,p) = (e,/i) or (e, r) or (fj,,r) . 



In general, these expressions are rather lengthy, but for 



•■/IT 



and J eM = —mi ^2 Amp, cf 2 sf 2 cf 3 sin 2a and I, 



- and U e3 = they simplify to 
— mi m 2 c\ ' 1 ! 



-12 s i2 s 23 sin 2a. 



The single physical low energy phase a corresponds to only one independent invariant. 
With the rephased mass matrix Eq. (fl4"]l we have = 0, J, 



-ml AB 2 D c 4 sin 20. By calculating the ratio I efM / I e 



and J eT : 

c 4 = cot 4 823- Moreover, we can write 



= -ml AB 2 D sin 20 
it follows again that 



mi m 2 Amg c 4 2 s\ 2 c 4 3 sin 2a ~ Am\ Am^ c 4 2 s\ 2 c 23 sin 2a 
= m\AB 2 D sin 20 = (m) m 2 , _B 2 D sin 20 , 



(16) 



which shows that is closely related to the low energy Majorana phase a. This phase a 
shows up in the effective mass governing neutrinoless double beta decay, which in our case 
is just [H] 

' ' ^ ~ (17) 



(m) 



1 — sin 2#i2 sin a 



We show in Fig. [T]a scatter plot of sin 2 9 12 against the effective mass for Eq. (El), where all 
parameters are complex. The oscillation parameters were required to lie in their current 3cr 
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ranges from Ref. [12J. One can identify basically two bands in this plot, one corresponding 
to (m) ~ a/ Am} cos 2^12 = a/ Am} (1 — 2 sin 2 #i 2 ), the other to (m) ~ a/ Am}. 
Let us compare our discussion with that in the /x-r symmetric case. With the generic 
prediction of 6*13 = 0, it follows that there is no Dirac phase. However, in contrast to the 
SSA case, in general both Majorana phases are present because there are in general three 
non-vanishing neutrino masses. This extra phase in /i-r case is hard to detect since as far 
as neutrinoless double beta decay is concerned, the effective mass will depend only on one 
phase due to the fact that U e3 = 0. The second phase is in principle observable in processes 
such as — and v T — v T oscillations p3] or rare decays such as K + — > n~ ji + [T3] , 
The experimental challenges in order to observe these processes at their predicted rates 
are however depressingly breathtaking, see the discussion in Ref. [T3] . 
For the case of L e — L^ — L T , there is also no Dirac phase, and in addition only one Majorana 
phase a, which is fixed by the symmetry to the value tt/2. The required (large and tuned) 
breaking of L e — L^ — L T will induce non-zero m 3 and #13 and therefore all three possible 
phases will be present. The phase a will stay close to 7r/2 and therefore the effective mass 
will stay close to (m) ~ a/ Am} cos2# 12 , unless the breaking parameters are as large as the 
parameters allowed by the symmetry. In case of SSA the effective mass can take any value 
between a/ Am} cos 2#i2 and a/ Am}, which can be used to distinguish the two scenarios. 
We stress however again that L e — L M — L T requires unnaturally large breaking. 

2.2 Cases (B) and (C) 

Let us now discuss the cases (B) and (C) given in Eq. ([3]) and Eq. (J3|). We note here that 
the effect of renormalization in case (B) is basically identical to case (A) discussed in the 
previous Subsection: c gets rescaled with (1 + 6 T ). In case (C) there is basically no running, 
because c gets rescaled with (1 + S^) <C (1 + 5 T ). In case (B) the mass matrix has the form: 
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Exactly analogous to case (A), in this case also we have a zero eigenvalue with an eigen- 
vector given by (— 1/a/1 + |c| 2 ,0,c/a/1 + |c| 2 ) t , therefore U^i = or = 0, depending 
on whether the hierarchy is normal or inverted. This is however incompatible with obser- 
vations which show that U^i 7^ for all i = 1,2,3. Similar situation holds for case (C) 
also, which predicts that U T \ = or U T s = 0. Therefore both these cases (B) and (C) are 
not phenomenologically viable. A possible way out would be to take correlations from the 
charged lepton sector into account, which however requires additional input and leads to 
less predictivity. 
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2.3 Breaking of Strong Scaling 



What happens if we break the SSA conditions? There are three possibilities for this, 
namely 
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As a consequence, non-zero and U e % will be generated and the prediction tan 2 2 3 = 1/ c 2 
receives corrections of order e. It turns out that in cases (Al) and (A2) \U e ^\ is roughly 
eD/2 and m3 is roughly De a/ Am 2 A /2. In case (A3) one finds jC^I ~ eD and m 3 is of 
order e 2 ^/ Am\. By allowing e to vary between zero and 0.22, we display the behavior of 
m 3 and \U e3 \ as a function of e in Fig. The oscillation parameters were again required 
to lie in their current 3a ranges from Ref. [12J. Only for sizable breaking of SSA of order 
of the Cabibbo angle one can probe such \Ues\ values in next generation experiments. The 
corrections to the effective mass given in Eq. ( TPTj) are of order 7713 \ U e s \ 2 and therefore 
completely negligible. 

2.4 Weak Scaling? 

One could relax the strong scaling condition and introduce a "weak scaling Ansatz" . This 
weaker version of SSA corresponds to the case where the equality holds for only a subset 
of flavors, e.g., 

These relations are also invariant under the usual renormalization group extrapolations. 
Note that breaking of SSA, as discussed in the previous Subsection, with large e corresponds 
to weak scaling. The third weak scaling condition has sometimes been invoked as a way to 
understand large atmospheric mixing simultaneously with small Am 2 & / Am\ [15] . Without 
additional input, for instance extra symmetries to make certain elements of A4 U equal, there 
is not much predictivity for these cases. Therefore, we do not consider this possibility 
anymore. 



3 Possible Theoretical Origin of Strong Scaling 

In this Section, we speculate on the theoretical origin of the strong scaling rule from fun- 
damental principles. We present two gauge models based on discrete family symmetries, 
which when broken lead to the strong scaling rule. 
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Scenario I: A gauge model 

First we present a gauge model, where a symmetry leads to the strong scaling rule. The 
model is based on the Standard Model gauge group supplemented by a family symmetry 
group -D4 x Z 2 as described in Ref. [9], where also the mathematical details of -D4 can be 
found. The leptons are assigned to the following representations of the discrete symmetry 
group: 

{L e , ( L l t )}~ ttf. 2+ ) ; {** ( ^ ) } ~ Or, 2-} ; (21) 

{(iv,),^),^)}-^-,^,^-}, 

where the subscripts +, — refer to the transformation under Z 2 and the rest are the D4 
representations. One requires [9] five Higgs doublets assigned to the group representations: 

0! ~ lr ; 4> 2 ~ it ; h ~ II ; ( ^ J ~ 2+ . (22) 

As has been shown in [9], this leads to the neutrino Dirac mass matrix of the form 

/ a \ 
m D = b d e 

\ / 



(23) 



a diagonal charged lepton mass matrix with arbitrary elements and also diagonal right- 
handed neutrino mass matrices with arbitrary elements. Using these, we can calculate the 
neutrino mass matrix using the type I seesaw formula to get 

/ a 2 A + b 2 B bdB beB \ 
M v = bdB d 2 B deB , (24) 

\ beB deB e 2 B J 

where A, B are the inverse values of the first two right-handed neutrino masses. Again, 
strong scaling with c = d/e is obeyed in this case. 

Scenario II: SSA from spontaneous breaking of \l—t symmetry 
Now we present another type I seesaw model in which SSA follows from spontaneous 
breaking of discrete /i-r symmetry embedded into a -D4 x Z 2 model. Consider an extension 
of the leptonic sector of the Standard Model with additional Higgs doublets H e fM T and H[ 2 
with transformation properties under D4 x Z2 given in Table [TJ The Yukawa couplings in 
this model can then be written as: 

£ = L e H e {hi N e + h 2 N fl + h 3 N T ) + (L M + L T H T ) {h 4 N e + h 5 N fM + h 6 N T ) 

+fn L e H e e R + /21 (Lfi Hf, + L T H T ) e R + / 22 {L^ fi R + L T tr) H[ ^5) 

+/ 33 (4 - L T t r ) H' 2 + \ £ M id NT C- 1 Nj + h.c. 
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Table 1: Transformation properties of the particle content of the model in Scenario IV. 



Symmetry breaking in general would lead to vevs of the form (H e ^ T ) = v e ^ T , all being 
different from each other. Using the seesaw formula, we can then calculate the light 
neutrino mass matrix }A V and find that it satisfies SSA with c = — . Note that scaling is 
obeyed independent of the form of M R . The atmospheric neutrino observations imply that 

~ v T which is a constraint on the model. 
Turning to the charged lepton masses, with the definition ip L Mg ijj R we can write 

( fllV e 

M e = /21 Vp f'22 v[ + /sa v' 2 

V J'2lV T 

where v[ 2 are the vevs of the H' 12 fields. For f 21 = 0, this can lead to the desired diagonal 
charged lepton masses with no contribution to U . For non-zero / 2 i one would expect non- 
zero U e 3 of the order of the 12 mixing angle in the matrix diagonalizing M^. If the entries 
of Mi are hierarchical, (M e ) 33 ~ m T 3> (M e ) 22 ~ m M , then 9[ 2 ~ (Mg) 2 i/m ll . 

4 Summary 

We have introduced a new concept of "strong scaling" to the neutrino mass matrix in 
which the ratio of certain elements of the neutrino mass matrix is a constant. We find it 
to be a new possibility to generate zero U e s and also the inverted hierarchy. In contrast 
to the flavor symmetry L e — — L T one can get correct phenomenology without any 
breaking of the Ansatz. Renormalization group running from the scale of neutrino mass 
generation to low scale has basically no effect on scaling, in particular U e3 = and m 3 = 
remain true for all energy scales. The phenomenology is predictive and testable, and the 
number of parameters equals the number of observables so that the Ansatz can be fully 
reconstructed. Observation of a departure of the atmospheric mixing angle from maximal 






/22 v[ - faa v' 2 



(26) 
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with tiny U e s would distinguish this model from an approximate /x-r symmetric mass 
matrix. Observation of a normal hierarchy for neutrinos or quasi-degenerate masses will 
rule out our Ansatz. We have also presented scenarios for the theoretical realization of the 
strong scaling Ansatz. 

Acknowledgments 

The work of R.N.M. was supported by the National Science Foundation grant no. Phy- 
0354401 and the Alexander von Humboldt Foundation (the Humboldt Research Award). 
The work of W.R. was supported by the "Deutsche Forschungsgemeinschaft" in the "Son- 
derforschungsbereich 375 fiir Astroteilchenphysik" and under project number RO-2516/3- 
1. R.N.M. would like to thank M. Lindner for discussions and hospitality at TUM during 
the time this work was done. We thank M. Ratz for comments. 

Note added: After the paper was submitted to the arXiv, it was brought to our attention 
that there is a model in the literature [16] leading to a mass matrix in our scaling Ansatz 
(A). It was also noted in Ref. [H] that the results m 3 = and U e3 = are invariant under 
RGE corrections. 

References 

[1] R. N. Mohapatr a et al, |hep-ph/0510213t R. N. Mohapatr a and A. Y. Smirnov, 
|hep-ph /0603 1 1 8[ A. Strumia and F. Vissani, [hep-ph /0606054, 

[2] P. Minkowski, Phys. Lett. B 67, 421 (1977); M. Gell-Mann, P. Ramond, and R. Slan- 
sky, Supergravity (P. van Nieuwenhuizen et al. eds.), North Holland, Amsterdam, 1980, 
p. 315; T. Yanagida, in Proceedings of the Workshop on the Unified Theory and the 
Baryon Number in the Universe (O. Sawada and A. Sugamoto, eds.), KEK, Tsukuba, 
Japan, 1979, p. 95; S. L. Glashow, The future of elementary particle physics, in Pro- 
ceedings of the 1979 Cargese Summer Institute on Quarks and Leptons (M. Levy et al. 
eds.), Plenum Press, New York, 1980, pp. 687; R. N. Mohapatra and G. Senjanovic, 
Phys. Rev. Lett. 44, 912 (1980). 

[3] S. T. Petcov, Phys. Lett. B 110 (1982) 245; an incomplete list of more recent studies 
is: R. Barbieri, L. J. Hall, D. R. Smith, A. Strumia and N. Weiner, JHEP 9812, 
017 (1998); A. S. Joshipura and S. D. Rindani, Eur. Phys. J. C 14, 85 (2000); 
R. N. Mohapatra, A. Perez-Lorenzana and C. A. de Sousa Pires, Phys. Lett. B 474, 
355 (2000); Q. Shafi and Z. Tavartkiladze, Phys. Lett. B 482, 145 (2000); K. S. Babu 
and R. N. Mohapatra, Phys. Lett. B 532, 77 (2002); H. J. He, D. A. Dicus and 
J. N. Ng, Phys. Lett. B 536, 83 (2002); G. K. Leontaris, J. Rizos and A. Psallidas, 
Phys. Lett. B 597, 182 (2004); P. H. Frampton and R. N. Mohapatra, JHEP 0501, 025 
(2005); S. T. Petcov and W. Rodejohann, Phys. Rev. D 71, 073002 (2005); W. Grimus 
and L. Lavoura, J. Phys. G 31, 683 (2005); G. Altarelli and R. Franceschini, JHEP 
0603, 047 (2006). 



11 



[4] T. Fukuyama and H. Nishiura, hep-ph/9702253; R. N. Mohapatra and S. Nussinov, 



Phys. Rev. D 60, 013002 (1999); E. Ma and M. Raidal, Phys. Rev. Lett. 87, 011802 
(2001); C. S. Lam, Phys. Lett. B 507, 214 (2001); PF. Harrison and W. G. Scott, 
Phys. Lett. B 547, 219 (2002); T. Kitabayashi and M. Yasue, Phys. Rev. D 67, 015006 
(2003); W. Grimus and L. Lavoura, Phys. Lett. B 572, 189 (2003) ; J. Phys. G 30 
73 (2004); Y. Koide, Phys. Rev. D 69, 093001 (2004); A. Ghosal, |hep-ph/0304090] 
W. Grimus, A. S. Joshipura, S. Kaneko, L .Lavoura, H. Sawanaka, M. Tanimoto, 
Nucl. Phys. B 713, 151 (2005); R. N. Mohapatra, SLAC Summer Inst, lecture; 
|http : //www-conf . slac . Stanford . edu/ssi/2004| JHEP 0410, 027 (2004); A. de 
Gouvea, Phys. Rev. D 69, 093007 (2004); S. Choubey and W. Rodejohann, Eur. 
Phys. J. C 40, 259 (2005); R. N. Mohapatra and W. Rodejohann, Phys. Rev. D 
72, 053001 (2005); R. N. Mohapatra and S. Nasri, Phys. Rev. D 71, 033001 (2005); 
R. N. Mohapatra, S. Nasri and H. B. Yu, Phys. Lett. B 615, 231 (2005); Phys. Rev. 
D 72, 033007 (2005); Y. H. Ahn et al. } Phys. Rev. D 73, 093005 (20 06); T. Ota and 
W. Rodejohann, Phys. Lett. B 639, 322 (2006); K. Fuki, M. Yasue, |hep^ph/0608042[ 

[5] R. N. Mohapatra in Ref. [I]. 

[6] S. T. Petcov and W. Rodejohann in Ref. [3]. 

[7] K. S. Babu, C. N. Leung and J. T. Pantaleone, Phys. Lett. B 319, 191 (1993); 
P. H. Chankowski and Z. Pluciennik, Phys. Lett. B 316, 312 (1993); S. Antusch, 
M. Drees, J. Kersten, M. Lindner and M. Ratz, Phys. Lett. B 519, 238 (2001). 

[8] G. C. Branco, R. Gonzalez Felipe, F. R. Joaquim, I. Masina, M. N. Rebelo and 
C. A. Savoy, Phys. Rev. D 67, 073025 (2003). 

[9] C. Hagedorn and W. Rodejohann, JHEP 0507, 034 (2005). 



[10] U. Sarkar and S. K. Singh, hep-ph/0608030 



[11] For a review, see C. Aalseth et al, hep-ph/0412300 



[12] T. Schwetz, [hep-ph/0606060] M. Maltoni et al, |hep-ph/0405172| y5. 



[13] A. de Gouvea, B. Kayser and R. N. Mohapatra, Phys. Rev. D 67, 053004 (2003). 

[14] W. Rodejohann, Ph ys. Rev. D 62, 013011 (2000); J. Phys. G 28, 1477 (2002); K. Zu- 
ber, |hep-ph/0008080 ; A. Atre, V. Barger and T. Han, Phys. Rev. D 71, 113014 (2005). 

[15] G. Altarelli and F. Feruglio, Phys. Lett. B 451, 388 (1999); for right-handed neutrino 
masses, mass matrices similar to those dictated by weak scaling have been considered 
by E. K. Akhmedov, G. C. Branco and M. N. Rebelo, Phys. Lett. B 478, 215 (2000). 

[16] W. Grimus and L. Lavoura, J. Phys. G 31, 683 (2005); Phys. Rev. D 62, 093012 
(2000); L. Lavoura, Phys. Rev. D 62, 093011 (2000). 



12 



0.055 




5- -035 f».\. V . * .• • V'- ... " 
.03 1 •-•.!* • v.*. * . ... ... r >.;.;... 




24 0.26 0.28 0.3 .0,32 0.34 0.36 0.38 0.4 

sin 2 12 



Figure 1: Scatter plot of sin 2 9 12 against the effective mass for Eq. 




Figure 2: Scatter plots of \U e3 \ and m 3 against e for broken SSA from Eq. ([191) . 



